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The XMM-Newton long look of NGC 1365: uncovering of 
the obscured X-ray source 
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ABSTRACT 

We present an analysis of the extreme obscuration variability observed during an 
XMM-Newton 5-days continuous monitoring of the AGN in NGC 1365. The source 
was in a reflection-dominated state in the first ^1.5 days, then a strong increase of the 
7-10 keV emission was observed in ~10 hours, followed by a symmetric decrease. The 
spectral analysis of the different states clearly shows that this variation is due to an 
uncovering of the X-ray source. From this observation we estimate a size of the X-ray 
source Ds < 10^^ cm, a distance of the obscuring clouds R ^ 10^^ cm and a density 
n ~ 10"'^^ cm~'^. These values suggest that the X-ray absorption/reflection originate 
from the broad line region clouds. This is also supported by the resolved width of the 
iron narrow Ka emission line, consistent with the width of the broad Hf3 line. 
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1 INTRODUCTION 

The structure, size and physical properties of the circumnu- 
clear medium of AGNs is still a matter of debate and inves- 
tigation. The axisymmetric (torus-like) geometry required 
by AGN unified models (e.g. Antonucci 1993) can be ob- 
tained by means of quite different structures, from galactic 
dust lanes to sub-parsec scale, high density tori. In this con- 
text, the hard X-ray observations of the nearby obscured 
Seyfert Galaxy NGC 1365 have proven to be a unique op- 
portunity to investigate the properties of the AGN circum- 
nuclear medium, as well as to put tight constraints on the 
size of the X-ray source. 

NGC 1365 is a nearby (z=0.0055) Seyfert 1.8 Galaxy 
Several observations performed in the hard X-rays in the 
past ~ 15 years with all the major X-ray observatories re- 
vealed several spectral changes from Compton-thin to reflec- 
tion dominated states, implying an extraordinary absorp- 
tion variability from Nh ~ 10^^ cm"^ to Nh > 10^* cm"^ 
(Risaliti et al. 2005A, 2007, hereafter R05A and R07). The 
most extreme variations have been observed during a recent 
Chandra campaign consisting of six short (15 ks) observa- 
tions performed in ten days. The source was in a Compton- 
thin state during the first observations, then "disappeared" 



(i.e. only the reflected component was visible) in the sec- 
ond observation two days later, and was back in the initial 
state during the third and later observations. The interpre- 
tation of this extreme variability is based on the presence of 
Compton-thick eclipsing clouds crossing the line of sight to 
the central source. We discussed this scenario in R07 where 
we obtained stringent constraints on the X-ray source size 
(D < 10^* cm) and, under the hypothesis of Keplerian mo- 
tion of the eclipsing clouds, on the distance of the absorber 
from the central black hole {R ^ 10^" cm). 
The main limits of this analysis are the lack of continuous 
monitoring during the eclipse, and the relatively low S/N of 
the Chandra spectra, due to the low effective area at high 
energies (E>2 keV) and the short observing time. In order 
to solve these problems, we obtained a 5-day XMM-Newton 
almost-continuous look of NGC 1365. During this obser- 
vation the source was on average in a low flux, reflection- 
dominated state. However, a dramatic spectral change was 
observed after ~ 2 days, when a hard transient component 
appeared for ~ 60 ks. In this paper we present the analysis 
of this XMM-Newton long look of NGC 1365, and we focus 
on the implications of the observed absorption variability. 
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2 DATA REDUCTION AND ANALYSIS 

The observation was performed from May 30 to June 5, 2007. 
for three complete orbits (1384-1386). The reduction and 
calibration were performed with the SAS package, following 
standard procedures as recommended by the XMM-Newton 
Science Operation Center. 

We included all the time intervals in the timing and 
spectral analysis, except for the final ~ 15 ks of the second 
orbit, when high background flares saturated the EPIC de- 
tectors. The spectra and light curves were extracted from 
a circular region with a 25 arcsec radius. With this choice, 
the background level is always lower than 5% of the source 
spectrum at all times and energy intervals. For both the tim- 
ing and spectral analysis, we always used the complete PN 
and MOS (merged) set of data, which are fully consistent 
with each other. In the figures however we only show the 
PN data, for clarity. 

In Fig. 1 we show the low-energy (2-5 keV) and high- 
energy (7-10 keV) light curves obtained from the PN. A 
visual analysis of these curves reveals a strong variation 
during the first half of the second orbit, consisting of an 
increase of the hard fiux, followed by a symmetric decrease. 
The hard fiux in the second half of the second orbit and in 
the third orbit remains constant, at a level slightly higher 
than in the first orbit. At low energies the source remains 
constant throughout the whole observation. Considering the 
past observed fast absorption changes, the most likely inter- 
pretation of this variability is a transient decrease of the 
absorbing column density during the second orbit. 

In order to test this hypothesis and to perform a com- 
plete spectral analysis of our data, we extracted the spec- 
tra from six time intervals, corresponding to different spec- 
tral states, as shown in Fig. 1 (the fifth and sixth interval 
show no obvious variations, and are separated only because 
they belong to two different orbits). The spectral analysis 
was performed with the XSPEC 12.4 package (Arnaud et 
al. 1996). 

Here we are only interested in the variations observed in 
the high-energy (E>2 keV) emission. Therefore, we do not 
discuss in detail the lower energy (0.5-2 keV) emission, which 
has a diffuse origin, as shown by Chandra images (R05A). 
However, this component cannot be completely neglected 
for two reasons: 1) an analysis of the Chandra diffuse emis- 
sion shows a faint, but not negligible, tail at energies higher 
than 2 keV; 2) the fits of the high energy spectra should 
not include components whose extrapolations towards low 
energies over-predict the observed emission. In order to take 
into account these aspects, we extracted a spectrum of the 
diffuse emission from the total Chandra data, in an annulus 
centered on the nucleus, with external radius equal to that 
used in the XMM-Newton analysis, and inner radius of two 
arcsec, in order to remove the central emission. We then fit- 
ted this spectrum with a multi-component model providing 
a good analytical representation of the data. The analysis 
of the XMM-Newton spectra was performed adding to all 
models the fixed components found in the Chandra analy- 
sis, with a 5% maximum range in the cross-calibration. 

We first performed the analysis of the six intervals sep- 
arately. The spectrum of the first interval can be entirely 
reproduced by a continuum reflection (PEXRAV model in 
XSPEC, Magdziarz & Zdziarski 1995) plus a set of three 
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Figure 1. 2-5 keV (top) and 7-10 keV (bottom) light curves of 
the nuclear emission of NGC 1365 during the XMM-Newton long 
look. The horizontal segments show the extraction time intervals 
for the six spectra discussed in the text. 

emission lines, at energies between 6.4 and 7.0 keV. The 
other five spectra are successfully reproduced with the same 
model as above, plus an intrinsic continuum component, 
consisting of an absorbed power law. The photon index, col- 
umn density and flux of the intrinsic component were left 
free to vary for each spectrum. The slope of the incident 
continuum for the reffection component is assumed to be the 
same as that of the intrinsic emission. In all cases the pa- 
rameters of the reffection continuum are compatible with the 
values found in the first spectrum. The column density and 
photon index of the intrinsic component are also constant 
within the errors, with values in the ranges F ~ 2.3 — 2.9 
and Nh ~ 5 — 10 x 10^^ cm^^. Instead, the fiux shows strong 
variations. 

As a second step in order to obtain a physically con- 
sistent view of the observed variability, we performed a si- 
multaneous analysis of the six spectra, imposing constant 
reffection parameters. The photon index and the absorption 
column density of the intrinsic power law were also required 
to be constant, while the flux was left free to vary in all the 
six intervals. The best flt parameters are shown in Tab. 1. 
In Fig. 2 we show the three spectra from the first, third and 
sixth interval, and the residuals for all the intervals. It is 
clear from the overall goodness of fit and from the residual 
distribution, showing no obvious unaccounted features, that 
the adopted model provides a good representation of the 
spectral variations observed in Fig. 1. 



3 DISCUSSION 

The spectral analysis presented above shows that the ob- 
served high energy variations are due to a power law com- 
ponent, absorbed by a column density Nh ~ 9 x 10'^'' cm~^, 
rising from ~zero up to its maximum fiux in ~ 30 ks, and 
then decreasing down to ~ 20% of its maximum in a sim- 
ilar time. These observations can be interpreted as due to 
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Table 1. XMM-Newton Analysis 



Energy (keV) 

Figure 2. Upper panel: PN Spectra obtained from the time in- 
tervals 1, 3 and 6 as shown in Fig. 1. The increase in the in- 
terval 3 spectrum (red diamonds) and, in a smaller fraction, in 
the interval 6 spectrum (green triangles) is well reproduced by 
an intrinsic power law continuum absorbed by a column density 
Nh ~ 9 X 10^3 cm-2 (see text for details). The other three PN 
and six MOS spectra are not shown for clarity, but are fully dis- 
cussed in the text. Middle panel: contribution to each of 
the six spectra, showing no significant feature not included in 
the model. Lower panel: best fit model components for the first, 
second and third interval. The continuous lines are common to 
all intervals, and represent the reflected components. The dashed 
lines show the intrinsic component, negligible in the first interval, 
showing up in the second interval, and further increasing in the 
third. 



intrinsic continuum variability, or to a change in the absorb- 
ing column density along the line of sight. The former case 
would imply a switch-on of the X-ray source from an almost 
completely reflection-dominated state (Tab. 1) in ~ 30 ks, 
followed by a strong decrease in a similar time. Such short 
time scales are unlikely in a disk-corona scenario, as dis- 
cussed in ROSA, however they cannot be completely ruled 
out. On the other side, absorption variability is a most likely 
scenario in this case, given the X-ray observational history 
of NGC 1365: 

- Unequivocal column density variations in short times scales 
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Table 2. Line Parameters 



Line 


E (keV) 


cr (eV) 




1 
2 
3 


6 64+° °2 

"•"^-0.04 
6 95+0 02 
°-^''-0.02 


OQ+14 

58-36 

< 91 


800 
180 
330 



(eV) 



+70 
-60 
+45 
-30 
+50 
-60 



" EWs are estimated with respect to the reflection component. 

(< 2 days) have been observed during the Chandra campaign 
(R07). Here we refer in particular to the Compton-thin vari- 
ations occurred during the last four observations, which can- 
not be interpreted as due to intrinsic continuum variability. 

- Another transit of an obscuring cloud has been detected 
in a previous XMM-Newton 60 ks observation (Risaliti et 
al. 2008, submitted). In this case the crossing cloud does 
not completely cover the source, and its column density is 
~ 3 X 102"^ cm~2. This makes the spectral analysis more 
complicated than in our case. However, the eclipsing time is 
of the same order as that observed here, while the underly- 
ing continuum remains constant. 

- Considering all the past X-ray observations where the 
source was in a Compton-thin state, strong variations of the 
primary continuum have never been observed: the estimated 
intrinsic flux is constant within a factor of at most two. 

Based on these arguments, we conclude that the in- 
terpretation invoking absorption variability is the most ap- 
propriate for our observations. Within this scenario, thick 
clouds move along the line of sight, uncovering and then 
covering again the primary X-ray source. The covering fac- 
tor of the thick clouds is almost complete during the flrst 
XMM-Newton orbit (INT 1 in Fig. 1), then decreases to 
a minimum, and then increases again, leaving only ~ 10% 
of the source uncovered in INT 5 and 6. It is not possi- 
ble from these observations to tell directly whether the un- 
covering is complete at the flux maximum. We note how- 
ever that the extrapolated 2-10 keV flux in INT 3 (F(2- 
10)~ 3 X 10"^^ erg s~^ cm~2^ adopting the continuum pa- 
rameters shown in Table 1) is compatible with the highest 
value measured in past observations (Risaliti et al. 2005B, 
hereafter R05B) . This is an indication that most of the X-ray 
source is probably uncovered during INT 3. 

Two schemes for the geometrical structure of the cir- 
cumnuclear absorber are possible (Fig. 3): 
1) The Compton-thick absorber is made of clouds with about 
the same size of the source. On average, few (1-3) of these 
clouds cross the line of sight. Depending on the fluctuations 
in the number and column density of the clouds the source 
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Figure 3. Possible models for the structure of the circumnuclear 
absorber. The black points represent Compton-thick clouds. The 
regions with a darker background in the "torus" are those where 
at least one thick cloud crosses the line of sight to the X-ray 
source (the red spot in the center). The regions with a lighter 
background show the Compton-thin lines of sight. 

can be seen in a reflection-dominated state for most of the 
time, with a short interval in a Compton-thin state, as ob- 
served. 

2) The thick clouds are bigger than the X-ray source by 
a factor of the order of the ratio between the continuous 
Compton-thick and Compton-thin intervals in our observa- 
tion. In this case our Hght curve can be reproduced by the 
crossing of a single thick cloud, followed by a short inter- 
val where only thin gas is present along the line of sight, 
followed by another "big" thick cloud. 

The observations discussed here can only constrain the 
properties of the thick clouds, but nothing can be directly 
inferred on the thin material. In principle, it can be in a dif- 
ferent physical state (either homogeneous, or distributed in 
clouds of different size) and distance (either closer or farther) 
from the central source. The variations of Compton-thin ab- 
sorption observed with Chandra (R07) in time scales of two 
days rule out a homogeneous material. Considerations on 
the ionization state (see below for more details) also rule 
out that the thin material is closer to the center than the 
thick clouds. However, a layer of thin clouds external to the 
thick ones, or a co-spatial mix of thick and thin clouds, are 
both possible. 

In order to estimate the relevant geometrical and phys- 
ical quantities (i.e. the size of the X-ray source, the size 
and distance from the center of the obscuring clouds, their 
density) we flrst note that the covering/uncovering time de- 
pends only on the source size and on the velocity of the 
thick clouds, but not on the dimensions of the obscuring 
clouds, nor on the structure of the thin gas. If the covering 
clouds (for the remainder of the analysis, we always refer to 
the Compton-thick component only) have different veloci- 
ties, the maximum size of the X-ray emitting region is still 
given by the velocity of the fastest cloud, v, times the eclips- 
ing time, T. In R07 we discussed possible ways to constrain 
these two quantities: 

- The time intervals between consecutive observations which 
caught the source in different states provided upper limits 
to the eclipsing time T <1 day. 

- Two main arguments provided constraints to the cloud 
velocity: 1) a limit to the ionization state of the absorbing 
gas, which cannot be too highly ionized in order to effec- 
tively absorb the primary component, and 2) a limit on the 



width of the narrow iron emission line, under the assumption 
that reflection and absorption are due to the same physical 
component. 

The new observations presented here allow an improved 
estimate of all the above constraints: 1) the eclipsing time T 
is directly measured from the observations (T ~10 hours); 2) 
the iron narrow emission line is now resolved and therefore 
the line width can be directly measured; 3) having a precise 
determination of the reflection-dominated spectrum and a 
measurement of the intrinsic continuum from past observa- 
tions (R05B) we can put a more stringent upper limit on 
the ionization parameter U of the obscuring cloud. Here we 
discuss the line width and ionization parameter arguments, 
following the scheme of R07. 

1. Ionization parameter. In order to obtain complete 
obscuration of the primary source, as observed in INT 1 
(Fig. 1), the obscuring gas cannot be too highly ionized. The 
intrinsic continuum can be estimated from previous X-ray 
observations. In particular, the XMM-Newton observation of 
January 2004, described in R05B, provided the lowest- A^i^, 
highest S/N spectrum available. We therefore adopted the 
best fit continuum model from this observation as the in- 
trinsic ionizing emission from the central X-ray source. We 
added this component to the best fit model for INT 1 and 
we allowed for ionized absorption, following the model by 
Done (1992). In order to reproduce the observed spectrum 
without over-predicting the continuum, the ionization pa- 
rameter must he U < 20 erg cm s~^, and the absorbing 
column density Nh > 2 x lO'^* cm^'^. From this limit we 
can obtain upper limits on the source size and on the cloud 
distance from the center, R, under the assumption of Kep- 
lerian cloud velocity v = {GMbh / R)^^^ , and assuming that 
the cloud size Dc is larger than the size of the X-ray source 
Ds- With simple algebra (see R07 for details) we obtain: 

R > {GMbhY^''[TLx /{Uma.NH)f''' (1) 

where Lx is the ionizing 1-150 keV radiation, which we take 
from R05B: Lx = 1.5 x lO''^ erg s"^ The central black hole 
mass, Mbh can been estimated through the relations be- 
tween the black hole mass and (1) the bulge velocity dis- 
persion [\og{MBH/MQ) = 7.3 ± 0.4(0.3) Oliva et al. 1995, 
Ferrarese et al. 2006], and (2) the bulge luminosity in the 
K-band [log(Msif/M0) = 7.8 ±0.4(0.3), Dong & De Rober- 
tis 2006, Marconi & Hunt 2003], where the two quoted errors 
refer to the statistical(systematic) effects. From these values 
we obtain R > 1.2 x 10^^ cm. The limits on the cloud ve- 
locity and the X-ray source size are v = (GMbh /R)^^^ < 
9,000M^^'' km/s, and Ds < 2.7 x lO^^M^''^ cm. This 
size can be expressed in units of gravitational radii, Rg ~ 
2GMbh/c^, obtaining Ds < 10 M~^'''Rg. 
2. Line width. The narrow emission lines present in all 
spectra are due to reflection by circumnuclear gas. It is 
reasonable to assume that the reflector and the absorbing 
clouds are the same physical component. This is supported 
by the occurrence of Compton-thick states in the past X- 
ray observations of NGC 1365, with a duty cycle of about 
1/3 (lyomoto et al. 2005, Risaliti et al. 2000, R05A, R07). 
It follows that the Compton-thick absorber covers ~l/3 of 
the solid angle, and therefore a significant reflection com- 
ponent must be produced by this gas. Interestingly, the ra- 
tio R between the normalizations of the refiection and in- 
trinsic continuum components (Table 1) in the uncovered 
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spectrum (INT 3) is about 1/3, further supporting the iden- 
tification between the reflector and the thick absorber. In 
this scenario, we can estimate the velocity of the obscuring 
cloud from the width of the narrow component of the neu- 
tral iron emission Hnc. Thanks to its large equivalent width 
in reflection dominated states, the long XMM-Newton ob- 
servation allows a precise analysis of the line profile, and an 
estimate of its width. The line is resolved at the 4 a level 
(FWHM(Fe)=3100llg„!J km/s. Table 2). The cloud velocity 
in the plane of the sky is then obtained by dc-projccting the 
measured line width in spherical symmetry. The final esti- 
mate of the cloud velocity v = 21001}3qq km/s. The linear 
dimension the X-ray source is then Ds ~ 6^4 x 10^^ cm. 
Using the black hole mass values discussed in the previ- 
ous paragraph, and assuming Keplerian motion of the ob- 
scuring cloud, we find - 2.0 ± 0.7 x lO^^M^/* cm and 
Ds = 2±l Mf 1 Rg. 

In order to have a physically reasonable size of the X-ray 
source, Ds ~ 10 Rq, the black hole mass must be of the or- 
der of 2-3 xlO® Mq. These values are consistent at a ~ 1 tr 
level with the mass estimate from the velocity dispersion- 
black hole mass correlation (see above in this Section). The 
estimate from the bulge luminosity-black hole mass relation 
is instead higher at a 2 cr level, but this is expected due 
to the strong starburst contribution to the bulge luminosity 
of NGC 1365. A further indication of a black hole mass of 
the order of a few l(f Mq comes from the Mbh-HP width 
correlation (Kaspi et al. 2005). We obtained the width of 
the Hl3 broad component from Schulz et al. 1999, while the 
optical luminosity of the AGN component has been esti- 
mated from the X-ray 2 keV intrinsic luminosity (R05B), 
assuming a standard AGN SED in the optical/UV (Risal- 
iti & Elvis 2004) and the optical/X-ray ratio of Steffen et 
al. (2006). With these values Mbh = 2 x 10^ Mq. This es- 
timate is rather uncertain, both on the assumptions on the 
SED and on the use of the Kaspi et al. correlation down to 
a luminosity lower than those for which the correlation has 
been tested. 



The above analysis, mainly focused on the size of the 
X-ray source, also suggests an identification between the 
thick circumnuclear refiector/absorber and the broad line 
region. First, we note that the resolved velocity of the iron 
emission line is is consistent with the width of the broad 
Hl3 emission line (FWHM(7i"/3)~ 1900 km s'^, Schulz et 
al. 1999), analogously to what recently found in another 
Seyfert Galaxy, NGC 7213, by means of Chandra grating 
spectroscopy (Bianchi et al. 2008). Further support to this 
indentification comes from the estimate of the size Dc and 
density n of the obscuring clouds. As shown by the configu- 
rations in Fig. 3, Dc cannot be uniquely measured from our 
analysis. However, a cloud size much higher than the source 
size is unlikely from a geometrical point of view. This is 
also suggested by the relatively short duration (< 2 days) of 
the eclipse observed with Chandra (R07). Assuming Dc=l- 
3xDs, and a column density Nh ~ 2 x 10^* cm~^ we obtain 
n ~l-3xl0^^ cm~^, in agreement with the estimates for the 
broad line clouds. 



4 CONCLUSIONS 

Our analysis clearly demonstrates two main points: 

1) The size of the X-ray source is Ds < 10^^ cm, corre- 
sponding to a few gravitational radii. 

2) The obscuring clouds have a distance from the center and 

a density typical of the broad line region. The measured iron 
emission line width is consistent with the width of the broad 
Hp line. This further supports an identification between the 
X-ray reflector/absorber and the broad line clouds. 

Our estimates rely on a simplified scheme, with ho- 
mogeneous spherical clouds moving with Keplerian velocity 
around the central source. Possible effects not considered in 
our analysis are a non-Keplerian component of the cloud 
velocity (for example, due to an outflow), non-spherical ge- 
ometries, and possible Nh variations within the obscuring 
cloud. All these considerations make the actual rmmbcrs un- 
certain within a small factor, but do not change the two 
general conclusions listed above. 
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